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1. Introduction 

The world average of a, is currently determined primarily by high-precision determinations 
on the lattice. However, in the past few years there has been a renewed interest in the precision de- 
termination of as from non-strange hadronic T decays. One reason for this is the recent calculation 
of the ^{ccf) coefficient in the perturbative contribution to the Adler function [|l]]. The improved 
precision of the perturbative contribution emphasizes the importance of bringing other systematic 
effects, previously treated with varying degrees of completeness, under control. In fact, a number 
of competing analysis methods now obtain results which are not fully consistent with one another, 
even when applied to the same data. 

There are several theoretical issues related to this discrepancy between different determina- 
tions. There is a question of which resummation scheme, contour-improved perturbation theory 
(CIPT) or fixed-order perturbation theory (FOPT), is best used in the analysis. This issue has 
been the focus of many recent analyses, and we have chosen to present the results of each scheme 
separately. The present analysis focuses on theoretical questions related to the non-perturbative 
contributions which, due to the relatively low value of the T mass, are not entirely negligible. We 
present a new framework that both maintains consistency between the dimension at which the op- 
erator product expansion (OPE) is truncated and the choice of sum rule weights, and provides a 
quantitative description of duality violations (DVs). A complete account of the work presented 
here can be found in Ref. [Q]. 



2. Theory 

The analysis of hadronic T decay begins with the ratio 

j,^^ nr-^v.^^<ironsl ^^^^ 

Experimentally, it is possible to decompose the non-strange contribution to this branching ratio 
into vector (V) and axial- vector (A) components. Given the currents Jy = u{x)'f-d{x) and 7^ = 
u{x)y^ 'f d{x) , the two-point current correlation functions are defined by 

n(;; EE ij d'xe'^-^{0\T{J^^{x)jll{m (2.2) 

= (/^^-<?V^)nt^i(^')+^^^^<l(^2) ^ ^2.3) 

where the superscripts (1,0) label spin. Defining the spectral functions Py\{s) = (l/;r)Im nY'^(s) 
with s = the invariant squared-mass of the hadronic system, the non-strange contributions to 
Eq. (2.1) can also be expressed as the sq = m\ version of the weighted integral ||3|] 

r2 «n / .\2 



(2.4) 



Since the 7 = contributions are i&[{mu zLmj)'^] suppressed and thus numerically negligible apart 
from the 71 contribution to p^\s), the differential versions or Eq. ( p^ ) provide experimental de- 
terminations of Py^^\s). 
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Figure 1: FESR contour The discontinuity 
across the cut is 2/Imn(s). 



Figure 2: Contour relating DVs to the spectral 
functions between s„,i„ < sq <oa_ 



Cauchy's theorem, applied to the contour in Fig. |l], impUes that, for any and any analytic 
weight w{s), Ily^^\s) satisfies the finite energy sum rule (FESR) relation [Q] 



so 



dsw{s)pl,[^°\s) 



(2.5) 



2KI J\s\=so 

The LHS of Eq. (2.5) can be determined experimentally for < fn^- This then allows QCD pa 
rameters such as a, to be related to experimental data, provided the OPE representation 



3 +■ 4 



(2.6) 



can be reliably employed on the RHS.^This requires not only that ^> J^qcd^ but also caution 
in treating contributions in the vicinity of the positive real s-axis, where the OPE is expected to 



breakdown for values of s below rn\. 



Defining the correction to the OPE caused by DVs as ^v,a{^) = Tly^*'^ {s) — Hy^^ {s), we may 



\OPEt 



relate the experimental spectral functions to the OPE by 

r dswis)pi'^'\s) = -^J dswis) [nZHs)+^v.Ais)] 

Jo Z.711 J \s\—S{) 

dsw{s)p^^^{s) , 



2ni J\s\=so 



(2.7) 
(2.8) 



where Py\ = {\/7t)^v,A{^)- The second line of Eq. ( p^ ) comes from an application of Cauchy's 
theorem using the contour in Fig. ^, provided Py'^{s), as expected, vanish fast enough as 5 — > oo 
[^]. To suppress DVs, previous analyses have typically employed only doubly -pinched weights, 
and often restricted their attention to sq = m\. With these precautions, all previous analyses have 
taken Py^ {s) to be negligible. 

3. Strategy 

In actual analyses, it is necessary to truncate the non-perturbative OPE series, both for practical 
reasons and because the series is likely asymptotic. Through the residue theorem, a term of order 

'Here we have ignored numerically negligible ^(m^ ^) dimension 2 contributions. We will also ignore the logarith- 
mic scale dependence of the C,-. For a more detailed discussion, see Ref. P]. 
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s" in the polynomial w{s) will pick out the term in the OPE of order 1/^+^. Consistency then 
requires that, if the OPE is truncated at 1^(1/^'"+'), the maximum degree of the polynomial weights 
employed should be « OPE contributions for weights of degree n will in general depend on the 
n + l OPE parameters a^, C4, Ce, . . ., Cin+i- As the number of such linearly independent weights 
is also n + l, one can not perform true fits using FESRs involving only a single value while 
maintaining such consistency. Our analysis therefore involves the range Smm < < m^, with Smin 
chosen such that theory gives a good description of the data above it. 

We will also consider the presence of DVs in this analysis. While quite generally one would 
like to understand the effect of DVs, the necessity of working with < makes their inclusion 
more relevant since at lower values the residual DV effects are expected to be larger. Since 
no systematic theory of DVs is currently available from QCD, we employ a physically motivated 
model to examine their contributions. This ansatz, based on Regge theory, large-A^c, and analyticity 
[|],is given by igll 

Pvli') = Kv.Ae-'''^-^^m{av^ + s^v,A) . (3.1) 

Introducing the DV ansatz above adds four new fit parameters per channel on top of the OPE 
parameters considered in earlier studies. The data needed to fit this expanded parameter set is 



generated by evaluating the LHS of Eq. ( |2.8| ) for a range of with a set of weights w{s). 

Including the ansatz, Eq. ( |3.1| ), in this analysis allows us to perform fits with polynomial 
weights w{s) that are less than doubly-pinched. In the following section, we first present the results 
of fits to the unpinched weight w{s) = 1. After demonstrating the ability of our model to describe 
unpinched sum rules, we then present results of fits involving multiple weights, both pinched and 
unpinched. By showing, in the latter fits, our model's ability to describe sum rules with a diverse 
variety of weights, we aim to establish that our model for DVs provides an accurate description of 
physics missing from the OPE representation of the correlators. 



4. Fits 



Two experiments have made public their T decay spectral functions, ALEPH [10, 11] and 



OPAL []12|]. The 2005 analysis of ALEPH is more recent and is based on more statistics, so it 
would be the natural choice to apply our fitting scheme. Unfortunately, the publicly available 
ALEPH data omits correlations due to unfolding. Since a reanalysis of this data by the ALEPH 
collaboration is ongoing, the present analysis will focus on the 1998 OPAL data. 

We first present the results of standard fits to Eq. (^) with data generated in the interval 
1.5 GeV^ < So < m^, where the lower bound is chosen based on fit quality and stability. Fig. || 
shows the results of fits to OPAL V channel data using the weight w{s) = I, which results in values 
for the strong coupling of 



a,(m2) =0.307 ±0.019 
a,(m2) =0.322 ±0.026 



(FOPT) , 
(CIPT) . 



(4.1) 
(4.2) 



The figures demonstrate not only the expected presence of duality violations in the V channel but 



also the ability of the ansatz, Eq. (3.1), to accurately describe the data. Plots of CIPT results are 
nearly identical to the FOPT results shown here. 
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0.14 




So (GeV-) s (GeV-) 



Figure 3: The result of FOPT fits to V channel data with w{s) = 1. Shown are, on the left, the OPE(+DV) 
and spectral integrals and, on the right, the the OPE( +DV) and OPAL experimental Py^^\s). Red curves 
give the full OPE+DV results, while blue curves show only the OPE contributions. 




1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



So (GeV-) s (GeV-) 

Figure 4: The result of FOPT fits to V and A channel data with w{s) = 1. Shown are, on the left, axial 
OPE( +DV) and spectral integrals and, on the right, the OPE( +DV) and OPAL experimental p^^^\s). Red 
curves give the full OPE+DV results, while blue curves show only the OPE contributions. 



A simultaneous fit to OPAL V and A channel data over the same sq interval using w{s) = I 
yields values for a,{m^^) of 0.308(18) and 0.325(25) for FOPT and CIPT, respectively. These 



values are completely consistent with the V-only results of Eq. ( [4.2[ ). Fig. ^ again demonstrates not 
only the expected presence of DVs in the A channel but also the ability of our model to accurately 
describe the data. Fits to the A channel alone, however, are less stable than the V-only case. 

To examine the ability of the ansatz. to describe sum rules involving pinched weights we must 



refine our strategy. Fitting the parameters of Eq. ( |3.1| ) using only a single pinched weight leads 
to unstable results due to the pinched weight's suppression of DVs. To proceed, we fit to data 
generated by more than one weight, always including the unpinched weight w{s) = \ to fix the 
parameters of the DV model. The strong correlations between data generated using a range of 
sq and different w{s), however, lead to correlation matrices with eigenvalues which are zero at 
machine precision, preventing a standard construction from being employed. To deal with this 
problem, we fit our set of OPE and DV parameters p by minimizing the alternate fit quality 

Q"=L I {lefi^o) -4:\4;p)) (cw)^:.' {&\4)-ii:\4-^p)) > (4-3) 
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Figure 5: The results of FOPTfits to V channel data with wo{s) ~ 1, W2{s) = 1 — (s/sq)'^, andw3{s) = (1 — 
s/so)'^{l + 2s /so). The left and right plots show the OPE( +DV} and spectral integrals for W2{s) and wj,{s), 
respectively. Red curves give the full OPE+DV results, while blue curves show only the OPE contributions. 

where Ie7^ is the weighted spectral integral, ij-^'' is the weighted theory integral which depends on 
the set of parameters p, and C^*^"^ is the full covariance matrix for the weight w. 

This approach allows us to study the ability of the ansatz to describe data simultaneously for a 
range of pinched and unpinched weights. of course differs from the standard function, as the 
latter includes also correlations between spectral integrals involving different weights. Being forced 
to work with means we lose the statistical interpretation of the distribution around the minimum 
of our fit quality present for the function. Care must then be taken to properly incorporate the 
effects of such "cross-moment" correlations on the errors and covariances of the fitted parameters. 
This is done here by performing a linear fluctuation analysis as detailed in Ref. [^. 

Fitting the OPAL V channel data over the range 1.5 GeV^ < Sq < m^ using the fit quality ( p3| ) 
and the set of weights wq{s) = 1, W2{s) = 1 — {s/sqY, and wt,{s) = (1 —s/sqY{\ + 2s/sq) yields 
values for a,{ml) of 0.304(19) and 0.322(31) for FOPT and CIPT, respectively. Fig. | shows the 
results of this analysis for the W2 and ws cases (the corresponding wo results are indistinguishable 
from those shown in Fig. ^). We find that our DV ansatz is able to provide an accurate simulta- 
neous description of the experimental spectral integrals for all of the weights employed, whether 
unpinched, singly-pinched, or doubly-pinched. Though not shown here explicitly, we find that an 
accurate simultaneous description is not possible using OPE contributions alone. The results for 
a, from the combined wo, W2, and W3 fit are clearly in excellent agreement with those from the 
single-weight, unpinched wq fit shown in Eq. (p!2[). While we have chosen to present the result of a 
fit using the specific basis of weights described above, a variety of different weight sets have been 
examined with excellent consistency throughout. 

5. Conclusions 

We have presented an improved strategy for FESR analyses of hadronic T decay data. First, 
the analysis enforces consistency between the degree of the weight employed and the dimension 
at which the OPE is truncated. Second, the analysis allows for a quantitative investigation of the 
impact of residual duality violations not present in previous analyses. This is accomplished by 
employing a model, Eq. (pTT]), for the duality violation contributions to the vector and axial-vector 
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channel spectral functions. By examining the stability of fits relative to both the range of energies 
used in the sum rules and the choice of weights, we conclude that our model is able to accurately 
describe the physics of duality violations in the vector and axial-vector channels. 

Our most stable results come from fits to only the vector channel data. Regardless, all of the 
results obtained from combined vector and axial-vector fits are consistent with the vector channel 
only results. We take as our main result the standard X'^' ^^(5) = 1, vector channel output of 
Eq. Q, a,{ml) = 0.307(19) for FOPT and a,{ml) = 0.322(26) for CIPT. For reference, the 
values of as{ml) obtained by OPAL from the same data are 0.324(14) and 0.348(21) for FOPT 
and CIPT, respectively. Our results are seen to have central values shifted downward by slightly 
more than 1 

signa, somewhat larger errors, and a somewhat smaller FOPT-CIPT difference. The larger errors 
are not unexpected, given the necessity of fitting four extra parameters per channel. 
Running our values up to the Z mass yields 



The slight asymmetry caused by scaling the errors has been averaged over. 
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a,{Ml) =0.1 169 ± 0.0025 
a^(M|) =0.1 187 ±0.0032 



{MS, nf = 5, FOPT) 
(Ms, nf = 5, CIPT). 
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